####################################################################################
###
###   Immigration, Public Housing and Support for the French National Front
###   Gloria Gennaro
###
###   Paper Fig A29
###
####################################################################################


library(ggplot2)
library(rdrobust)

################################################################################
# Set Up
################################################################################

# Load Working Directories
source(paste0(wd_main, '/2_code/00_working_directories.R'))

# Load data and clean
source(paste0(wd_code, '/02_data_load_clean_rdd.R'))


data = df_did

################################################################################
# Bdw select
################################################################################

rdbwselect(data$hlm_share_log, data$running_pre, fuzzy=data$sruT, c=0) %>% summary()
df = data[abs(data$running_pre)<898,]

controls = "hlm_share_99 + res_owned_share_99 + imm_share_99 + fn_1995"

################################################################################
# Effect on immigration  (Figure A29a)
################################################################################

df['y'] = lm(as.formula(paste0('imm', " ~ ", controls)), data=df)$residuals
rdplot(df$y, df$running_pre, p=2, x.label='Population around the Cut-off', y.label='', title='')
ggsave(paste0(wd_res, '/figures/FigA29a.pdf'), height = 7, width = 5)


################################################################################
# Effect on home ownership  (Figure A29b)
################################################################################
df['y'] = lm(as.formula(paste0('res_own_share', " ~ ", controls)), data=df)$residuals
rdplot(df$y, df$running_pre, p=2, x.label='Population around the Cut-off', y.label='', title='')
ggsave(paste0(wd_res, '/figures/FigA29b.pdf'), height = 7, width = 5)
